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I. INTRODUCTION 


One of the most exciting recent developments in particle physics is the discovery of 
neutrino oscillations pointed out by the Super-Kamiokande experiment [1] and confirmed 
by the Sudbury Neutrino Observatory [2]. Neutrino oscillations provide evidence for non¬ 
vanishing neutrino masses and mixings, with the novel feature that large leptonic mixing 
angles are required, in contrast to what happens in the quark sector. Indeed, the combined 
results from these experiments suggest that, in addition to the large mixing angle required 
by the atmospheric neutrino data, another large angle should be present in the leptonic 
sector. This leads to the so-called large mixing angle (LMA) solution of the solar neutrino 
problem which turns out to be presently the most favored scenario for the explanation of 
the solar neutrino deficit. From a theoretical point of view, understanding the large leptonic 
mixing is still an unresolved mystery for which a considerable number of solutions have been 
proposed [3]. On the other hand, the appearance of neutrino masses much smaller than 
those of charged leptons is elegantly explained through the seesaw mechanism [4] which can 
be implemented by extending the standard model (SM) particle content with right-handed 
neutrinos. These can be easily accommodated in grand unified theories (GUT) where they 
appear on equal footing with the other SM particles. 

The heavy singlet neutrino states can also play an important role in cosmology, namely, in 
the explanation of the observed cosmological baryon asymmetry. During the last few years, 
the data collected from the acoustic peaks in the cosmic microwave background radiation [5] 
has allowed to obtain a precise measurement of the baryon asymmetry of the universe (BAU). 
The MAP experiment [6] and the PLANCK satellite [7] planned for the near future should 
further improve this result. At the present time, the measurement of the baryon-to-entropy 
ratio Y b = n B /s is 

0.7 x 1CT 10 < Y b < 1.0 x 1CT 10 . (1) 

Leptogenesis is one of the most attractive mechanisms to generate the BALL As first 
suggested by Fukugita and Yanagida [8], the key ingredient in leptogenesis are the heavy 
Majorana neutrinos which, once included in the SM, can give rise to a primordial lepton 
asymmetry through their out-of-equilibrium decays. This lepton asymmetry is subsequently 
reprocessed into a net baryon asymmetry by the anomalous sphaleron processes. 
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In spite of being attractive and successful, leptogenesis turns out to be extremely difficult 
or even impossible to test experimentally in a direct way. This difficulty is obviously related 
to the large masses of the heavy Majorana neutrino singlets. Nevertheless, the joint analysis 
of leptogenesis and low-energy neutrino phenomenology can be viewed as an indirect way 
of testing it and here the experimental results from neutrino oscillation experiments such 
as those related to the search of leptonic CP violation in the future long-baseline neutrino 
experiments are extremely valuable [9]. 

In this paper, we will address the question of linking the amount and sign of the BAU to 
low-energy neutrino experiments, namely to the sign and strength of the CP asymmetries 
measured through neutrino oscillations. Our analysis is performed in the weak basis (WB) 
where the charged lepton mass matrix rri( and the right-handed Majorana matrix Mr are 
both real and diagonal. In this WB, all CP-violating phases are contained in the Dirac 
neutrino mass matrix m/> The matrix me is arbitrary and complex, but since three of its 
nine phases can be eliminated through rephasing, one is left with six independent physical 
CP-violating phases. In order to study the link between the BAU generated through lepto¬ 
genesis and CP violation at low energies, it is crucial to use a convenient parametrization of 
rrifj. We shall make use of the fact that any arbitrary complex matrix can, without loss of 
generality, be written as the product of a unitary matrix U and a lower triangular matrix 
Ya. We show that U contains three phases which do not contribute to leptogenesis, while 
the other three phases contained in Ya contribute to both leptogenesis and low-energy CP 
violation. As a result, a necessary condition for having a link between leptogenesis and 
low-energy CP breaking is that the matrix U contains no phases, the simplest choice being 
obviously U — U. Within this class of Dirac neutrino mass matrices, we perform a search of 
the minimal scenarios where not only a good fit of low-energy neutrino data is obtained but 
also a link between the observed size and sign of the BAU and the strength of CP violation 
observable at low energies through neutrino oscillations can be established. 

II. GENERAL FRAMEWORK 

We work in the framework of a minimal extension of the SM which consists of adding to 
the standard spectrum one right-handed neutrino per generation. Before gauge symmetry 
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breaking, the leptonic couplings to the SM Higgs doublet 0 can be written as: 


Cy = ~Y„ (l», t>£) 4>v° r - Y, ( 4 °, i>») ^ 4 + H.c. 


( 2 ) 


where 0 = ?'t 20*. After spontaneous gauge symmetry breaking the leptonic mass terms are 
given by 


r = 

J—'m. 


Vl^dVr + -v° r t CM r v° + £°m e 


+ H.c. 


^n T L CM*n L + i°m e i 0 R 


+ H.c., 


( 3 ) 


where mo = vY u is the Dirac neutrino mass matrix with v = (0°)/\/i 2 ~ 174GeV, Mr and 
rri( = v Yp denote the right-handed Majorana neutrino and charged lepton mass matrices, 
respectively, and tir = (z/£, (v R ) c ). Among all the terms, only the right-handed neutrino 
Majorana mass term is SU(2) x U(l) invariant and, as a result, the typical scale of Mr can 
be much above the electroweak symmetry breaking scale v, thus leading to naturally small 
left-handed Majorana neutrino masses of the order m 2 D /MR through the seesaw mechanism. 
In terms of weak-basis eigenstates the leptonic charged current interactions are given by: 


£ w = -^wyp i »r4+H.c.. (4) 

It is clear from Eqs. (3) and (4) that it is possible to choose, without loss of generality, a weak 
basis (WB) where both nip and Mr are diagonal, real and positive. Note that in this WB, 
rrio is a general complex matrix which contains all the information on CP-violating phases. 
Since in the present framework there is no A L = 2 mass term of the form CMru° L) the 

total number of CP-violating phases for n generations is given by n(n — 1) [10] which are 
all contained in mo in this special weak basis 1 . 

We recall that the full 6x6 neutrino mass matrix A4 is diagonalized via the transforma¬ 
tion: 


V t M*V = V, (5) 

where V = diag(mi, m 2 , m 3 , Mi, M 2 , M 3 ), with 774 and M t denoting the physical masses of 
the light and heavy Majorana neutrinos, respectively. It is convenient to write V and T> in 

1 The counting of independent CP-violating phases for the general case, where besides mo and Mr there 
is also a left-handed Majorana mass term at tree level has been discussed in Ref. [11]. 


4 



the following form, together with the definition of M. : 


V = 




M 


( 0 m D 
rrio M r 


From Eq. (5) one obtains, to an excellent approximation, the seesaw formula: 


d v ~ -iC f rn D M n l ml K* = K ] M v K* , 


( 6 ) 


( 7 ) 


where M. v is the usual light neutrino effective mass matrix. The lcptonic charged-current 
interactions are given by 

~^= (4 7#1 Kv l + 4 Q N l ) + H.c., (8) 

where and N t denote the light and heavy neutrino mass eigenstates, respectively. The 
matrix K which contains all information on mixing and CP violation at low energies can 
then be parametrized, after eliminating the unphysical phases, by K = UsP with P = 
diag(l, e ta , e l/3 ) (i a and (3 are Majorana phases) and Us a unitary matrix which contains 
only one (Dirac-type) phase 5. In the limit where the heavy neutrinos exactly decouple 
from the theory, the matrix K is usually referred as the Pontecorvo-Maki-Nakagawa-Sakata 
mixing matrix, which from now on we shall denote as U u . 


CP violation in neutrino oscillations 


It has been shown [12] that the strength of CP violation at low energies, observable for 
example through neutrino oscillations, can be obtained from the following low-energy WB 
invariant: 

T CP = Tr [ H v , H e ] 3 = 6i A 21 A 32 A 31 Im [ (^)i 2 (^) 23 (^) 3 i ] , (9) 

where H u = ml and A 21 = (m M 2 — m e 2 ) with analogous expressions for 

A 31 , A 32 . This relation can be computed in any weak basis. The low-energy invariant (9) 
is sensitive to the Dirac-type phase 6 and vanishes for <5 = 0. On the other hand, it does 
not depend on the Majorana phases a and f3 appearing in the lcptonic mixing matrix. The 
quantity Tcp can be fully written in terms of physical observables once 

Im [ (7"4) 12 (7"4) 23 (7"C) 3 i ] = -A777,2! Am 32 J C p , (10) 
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where the Arne’s are the usual light neutrino mass squared differences and Jcp is the 
imaginary part of an invariant quartet appearing in the difference of the CP-conjugated 
neutrino oscillation probabilities P{v e —» /y /t ) — P(u e —> P M ). One can easily get: 


Jcp = Im [ {U v ) u {U v ) 22 {U v )\ 2 {U v y 21 } 

1 

= - sin(2 9 1 2 ) sin(2 d i3 ) sin(2 0 23 ) cos($i 3 ) sin <5 , 


( 11 ) 


where the 9ij are the mixing angles appearing in the standard parametrization adopted in 
[13]. Alternatively, one can use Eq. (10) and write: 

Im[(^) 12 (^) 23 (^) 31 ] 


Jcp = 


A m .21 Ami A m| 2 


( 12 ) 


This expression has the advantage of allowing the computation of the low-energy CP invari¬ 
ant without resorting to the mixing matrix U v . 

It is also possible to write WB invariants useful to lcptogenesis [12] as well as WB invariant 
conditions for CP conservation in the leptonic sector relevant in specific frameworks [11, 14], 

III. CP ASYMMETRIES IN HEAVY MAJORANA NEUTRINO DECAYS 


The starting point in lcptogenesis scenarios is the CP asymmetry generated through the 
interference between tree-level and one-loop heavy Majorana neutrino decay diagrams. In 
the simplest extension of the SM, such diagrams correspond to the decay of the Majorana 
neutrino into a lepton and a Higgs boson. Considering the decay of one heavy Majorana 
neutrino Nj, this asymmetry is given by: 

T(Nj -> £</>)- T(Nj -y£^) 


J = 


T(Nj -> t(f>) + T(Nj -> £(f) t) 

In terms of the Dirac neutrino Yukawa couplings the CP asymmetry (13) is [15]: 

1 


£ i = 


8 'ir(YvY I/ /jj 

where the index j is not summed over in (' YJY v )jj . The loop function f(x) includes the 
one-loop vertex and self-energy corrections to the heavy neutrino decay amplitudes, 


E Im[«U)i ]/ (M2 


M£ 


(13) 


(14) 


fix) = 


(1 + x) In 


x 


+ 


2 — x 


(15) 


1 + X J 1 — X 

From Eq. (14) it can be readily seen that the CP asymmetries are only sensitive to the 
CP-violating phases appearing in YjY u (or equivalently in rrv D m D ) in the WB where Mu 
and rri( are diagonal. 
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A. Hierarchical case: M\ < M2 Ms 


In the hierarchical case M\ < M 2 <C M 3 , only the decay of the lightest heavy neutrino 
N\ is relevant for lcptogenesis, provided the interactions of N { are in thermal equilibrium at 
the time 7V 2)3 decay, so that the asymmetries produced by the latter are erased before Ni 
decays. In this situation, it is sufficient to take into account the CP asymmetry gq. Since 
in the limit x 3> 1 the function f(x) can be approximated by 2 f(x) ~ — 3/(2y / a;), we have 
from Eq. (14) 


£1 = 


i6tt (yjyjn 23 

which can be recast in the form [16] 


£ M(buA ,Ml 




£ 3 Mi Im [ Yj Y u D^ 1 Yj Y* ] u 3M, Im [Yj M „ Y v * ] n 

£l “ 16 vr (yjyjn I 6 vru 2 (yjyjn 

using the seesaw relation given in Eq. (7). 


(16) 


(17) 


B. Two-fold quasi-degeneracy: Mi ~ M2 Af 3 

In the context of thermal lcptogenesis, when the observed baryon asymmetry is generated 
through the decays of the lightest heavy Majorana neutrino N\, there exists an upper bound 
on the CP asymmetry £1 which directly depends on the mass of the lightest neutrino M\ 
[17]. In turn, such a bound implies a lower bound on the lightest mass Mi, typically Mi > 
10 8 GeV. The latter bound 3 is however barely compatible with the reheating temperature 
bound T r{ < 10 8 ~ 10 9 GeV required in several supergravity models in order to avoid a 
gravitino overproduction [19]. To overcome this problem 4 * * * one can consider, for instance, the 
decays of two heavy neutrinos which are quasi-degenerate in mass, Mi ~ M 2 . In this case, 
the CP asymmetries £ t are enhanced due to self-energy contributions [ 21 ] and the required 
baryon asymmetry can be produced by right-handed heavy neutrinos with moderate masses 

2 This approximation can be reasonably used for x > 15. 

3 A more stringent constraint, Mi > 10 10 GeV, is obtained in [18]. 

4 Another possible solution to the gravitino problem is to consider non-thermal production mechanisms 

[20]. Since in these cases the condition Mi < Tr is not required, the gravitino problem is easily avoided 

once heavy particles can be created with a relatively low reheating temperature without threatening big 

bang nucleosynthesis. 
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M\ ~ M 2 < 10 8 GeV. Moreover, it has been shown that in the presence of small Dirac-type 
leptonic mixing at high energies and GUT-inspired Dirac neutrino Yukawa couplings, the 
heavy Majorana neutrino degeneracy is compatible with the LMA solar solution [22], 

Let us assume that the heavy Majorana neutrinos N\ and 1V 2 are quasi-degenerate. It 
is useful to define the parameter 5n which represents the degree of degeneracy between the 
masses M\ and M 2 as 


A 1 


(18) 


Since Mi ~ M 2 , we expect Sn -C 1. For the perturbative approach to remain valid, the 
tree-level decay width r, for each of the heavy Majorana neutrinos must be much smaller 
than the mass difference between them. This is translated into the relations 


D = [ Hv l U Ml « M 2 - Ml = 5 n Ml , i = 1,2, 

8 71 

where H u = YJ Y u . From this equation we can find the following lower bound for Sn- 

(■ H u )a Mi 


5 m Y> max 


(19) 


( 20 ) 


i=l,2 


8 7r Mi 

Assuming that this criterium is verihed, the CP asymmetries £i can be obtained combining 
Eqs. (14) and (18). We find 


* - { Im 1 (H ^ mi ] g, 

^ = - 57(^{ to[(H ' ) “ 1/[(1 + fw ) " 2 l - i Im[(Jf ' ) » 1 S - (21) 

Taking into account that for <5 n -C 1 the function /[ (1 + (bv) ±2 ] can be approximated by 5 : 

/[(1+<M ±2 ]~W-, (22) 

2 On 

we obtain: 


£ j = 


Im [ (H, 


2 ■ 

V) 21 


+ 3 Im l 1 jt 


7 = 1,2 


(23) 


16 7T (H v )jj L 5 N - 1V1 3 _ 

Typically, the term proportional to Mj/M 3 can be neglected and in this case £1 and £2 have 
the same sign. This aspect turns out to be relevant for the discussion on the relative sign 
between the BAU and low-energy leptonic CP violation. 


5 For Sn < 10 2 the error associated to this approximation is less than 1%. 



IV. ON THE CONNECTION BETWEEN LEPTOGENESIS AND LOW-ENERGY 
CP VIOLATION 

In this section we analyze the possible link between CP violation at low energies, mea¬ 
surable for example through neutrino oscillations, and leptogenesis. The possibility of such 
a connection has been previously analyzed in the literature [12, 23, 24], Nevertheless, we 
find it worthwhile presenting here a thorough discussion on the subject. In particular, we 
will address the following questions: 

• If the strength of CP violation at low energies in neutrino oscillations is measured, 
what can one infer about the viability or non-viability of leptogenesis? 

• From the sign of the BAU, can one predict the sign of the CP asymmetries at low 
energies, namely the sign of Jcp ? 

We will show that having an explicit parametrization of irtr, (or equivalently of Y v = mp/v) 
is crucial not only to determine which phases are responsible for leptogenesis and which ones 
are relevant for lcptonic CP violation at low energies, but also to analyze the relationship 
between these two phenomena. 

From the available neutrino oscillation data, one obtains some information on the effective 
neutrino mass matrix M. v which can be decomposed in the following way: 

U v d v Uj = M u = LL t , L = i m D D~ 1/2 . (24) 

The extraction of L from JA V suffers from an intrinsic ambiguity [25] in the sense that, given 
a particular solution Lq of Eq. (24), the matrix L = LqR will also satisfy this equation, 
provided that R is an arbitrary orthogonal complex matrix, R e 0(3, C), i.e. R R 7 =1. It 
is useful to take as a reference solution L 0 = U v dJ , so that: 

L = U v d]/ 2 R . (25) 

Since three of the phases of rrifj can be eliminated, the matrix L has 15 independent pa¬ 
rameters. The parametrization of L given in Eq. (25) has the interesting feature that all its 
parameters are conveniently distributed among U v , d v and R , which contain 6 (3 angles + 
3 phases), 3 and 6 (3 angles + 3 phases) independent parameters, respectively. Of the 18 
parameters present in the Lagrangian of the fundamental theory described by rrio and Dr, 
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only 9 appear at low energy in M. v through the seesaw mechanism. To further disentangle 
mo from Dr in L, one needs the 3 remaining inputs, namely the three heavy Majorana 
masses of Dr. As for the meaning of the information encoded in R , it turns out that the 
pattern of this matrix has a suggestive interpretation in terms of the different roles played 
by the heavy neutrinos in the seesaw mechanism. In fact, R can be viewed as a dominance 
matrix [26] since it gives the weights of each heavy Majorana neutrino in the determination 
of the different light neutrino masses rn t [27]. The fact that R^ are weights for m, is quite 
obvious clue to the orthogonality of R: 


rrii = 




(26) 


On the other hand, since Uj mo = —i d\! 2 RD l J . 2 , the single contribution rn t R,f: is also given 




by: 


■miR^ = 


{Vim 


D 


M, 


r Y 
'ji — yv tj 

~ ' - M, 


(27) 


Therefore, once U v is fixed, each weight Kj- just depends on the mass Mj of the j-th heavy 
Majorana neutrino and on its couplings with the left-handed neutrinos {mo)kj- Thus, the 
contribution of each heavy neutrino to rrii is well defined and expressed by the weight Re(/?7). 
One may roughly say that the heavy Majorana neutrino with mass Mj dominates 6 in m, if 


|Re(X y )| >> \Re(X ik )\ 


M~ 


Mi, 


k ^ j 


(28) 


which implies |Re(i??)| S> |Re(i?^)| . So that, if one of the heavy Majorana neutrino 
neutrino gives the dominant contribution to one of the masses rrii, this information is encoded 
in the structure of R. The interpretation in terms of weights is straightforward for the 
rotational part of R. However, one has to be careful because in the presence of the three 
boosts (controlled by the three phases) the weights R'k are not necessarily real and positive. 
Although this situation is more subtle, the above dominance arguments still hold. 

Coming back to the connection between lcptogenesis and low-energy data, it is important 
to note that U v does not appear in the relevant combination for leptogenesis YJY V , in the 
same way as R does not appear in M. u . Indeed, one has: 


m' D m D 


= D){ 2 R ] d u RD 


1/2 

R 


(29) 


6 See for instance Refs. [29, 30] for other approaches to the dominance mechanism. 
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From the above discussion, it follows that it is possible to write mo in the form mo — 
—i U u d l J 2 R D l J, 2 in such a way that leptogenesis and the low-energy neutrino data (contained 
in M. v ) depend on two independent sets of CP-violating phases, respectively those in R and 
those in U v . In particular, one may have viable leptogenesis even in the limit where there 
are no CP-violating phases (neither Dirac nor Majorana) in U v and hence, no CP violation 
at low energies [28]. Therefore, in general it is not possible to establish a link between 
low-energy CP violation and leptogenesis. This connection is model dependent: it can be 
drawn only by specifying a particular ansatz for the fundamental parameters of the seesaw, 
mo and Dr, as will be done in the following sections. 

The relevance of the matrix R for leptogenesis can be rendered even more explicit [27] 
by rewriting the e 1 asymmetry by means of Eq. (29) and defining R tJ = \R, l j\e v ' Pi:, l' 2 , A m 2 Q = 
A an( l Am| = Am| 2 - hi the case of hierarchical heavy Majorana neutrinos, say M\ <C 
M 2 -C M 3 one obtains 

3 Mi Am||_R 3 i | 2 sin <^ 3 i - Am%\R u \ 2 sm(p u 
£ i — —- Y - ~ 17? 9 i - 1 p | 2 T -, (30) 

lo 7 r v z mi|/inp + m 2 |/t 2 i| + m 3 |/£ 3 ip 


and we recover what one would have expected by intuition, namely that the physical quanti¬ 
ties involved in determining sq are just Mi, the spectrum of the light neutrinos, m t , and the 
hrst column of R, which expresses the composition of the lightest heavy Majorana neutrino 
in terms of the light neutrino masses m*. In the case Mi ~ M 2 -C M 3 , similar expressions 
hold, 

1 Mi Ad 2 lm[(R^d u R) 2 21 ] 

£j - 167 TV 2 M 2 - Ml {WdvR)^ 

where now also the mass M 2 and the second column of R are involved. A detailed study of 
the relevance of the matrix R for leptogenesis is under way [31]. 

As stressed before, different ansdtze for R have no direct impact on CP violation at low 
energy; the impact is in a sense indirect because R specifies if dominance of some heavy 
Majorana neutrino is at work in the seesaw mechanism [26]. 

In conclusion, the link between leptogenesis and low-energy CP violation can only be es¬ 
tablished in the framework of specific ansdtze for the leptonic mass terms of the Lagrangian. 
We shall derive a necessary condition for such a link to exist. In order to obtain this con¬ 
nection, it is convenient to use a triangular parametrization for mo, which we describe 
next. 
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Triangular parametrization 


It can be easily shown that any arbitrary complex matrix can be written as the product 
of a unitary matrix U with a lower triangular matrix Y A - In particular, the Dirac neutrino 
mass matrix can be written as: 


rn D = vUY a , 


(32) 


with Y a of the form: 


>A = 


2/21 e ^ 21 

v 2/31 e^ 31 


0 0 ^ 

2/22 o 

2/32 e 1 ^ 32 2/33 / 


(33) 


where r/ij are real positive numbers. Since U is unitary, in general it contains six phases. 
However, three of these phases can be rephased away by a simultaneous phase transformation 
on is®, P®, which leaves the lcptonic charged current invariant. Under this transformation, 
rrifj —> Ppmo , with P t t = diag (e*^ 1 , e*^ 2 , e*^ 3 ). Furthermore, Y A defined in Eq. (33) can be 
written as: 


U A = P% Y a Pp 


(34) 


where Pp = diag(l, e l/?1 , e 1 ^ 2 ) with 0 i = — 0 2 i, 02 = — 03 i and 


Y/\ = 


^ 2 /n 0 0 ^ 

2/21 2/22 0 

y 2/31 2/32 e ia 2/33 J 


(35) 


with a = 0 32 — 03i + 02 i- It follows then from Eqs. (32) and (34) that the matrix rri£> can 
be decomposed in the form 


m D = v Up P a Y a Pp , (36) 

where P a = diag(l, e* ai , e ia2 ) and U p is a unitary matrix containing only one phase p. 
Therefore, in the WB where rrq and Mr are diagonal and real, the phases p, aq, aq, ct, 
0i and 02 are the only physical phases characterizing CP violation in the leptonic sector. 
The phases relevant for leptogenesis are those contained in m' D m D . From Eqs. (34)-(36) we 
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conclude that these phases are cr, f3\ and (3 2 , which are linear combinations of the phases 
( f>ij . On the other hand, all the six phases of rrio contribute to the three phases of the 
effective neutrino mass matrix at low energies [12] which in turn controls CP violation in 
neutrino oscillations. Since the phases aq, a 2 and p do not contribute to leptogenesis, it is 
clear that a necessary condition for a direct link between leptogenesis and low-energy CP 
violation to exist is the requirement that the matrix U in Eq. (32) contains no CP-violating 
phases. Note that, although the above condition was derived in a specific WB and using the 
parametrization of Eq. (32), it can be applied to any model. This is due to the fact that 
starting from arbitrary leptonic mass matrices, one can always make WB transformations 
to render rrq and Mr diagonal, while rn R has the form of Eq. (32). A specific class of 
models which satisfy the above necessary condition in a trivial way are those for which 
t/ = l, leading to me = vY&. This condition is necessary but not sufficient to allow for 
a prediction of the sign of the CP asymmetry in neutrino oscillations, given the observed 
sign of the BAU together with the low-energy data. Therefore, we will consider next a more 
restrictive class of matrices rriu of this form and we will show that, in an appropriate limiting 
case, our structures for uid lead to the ones assumed by Frampton, Glashow and Yanagida 
in [24], 

V. MINIMAL SCENARIOS 

From the analysis carried out in the previous section, it becomes clear that the computa¬ 
tion of the cosmological baryon asymmetry Y R in leptogenesis scenarios strongly depends on 
the Yukawa structure of the Dirac neutrino mass term and on the heavy Majorana neutrino 
mass spectrum. Moreover, if one assumes that the seesaw mechanism is responsible for the 
smallness of the neutrino masses, then the connection between the baryon asymmetry and 
low-energy neutrino physics is unavoidable. In fact, this constitutes an advantage for the 
leptogenesis mechanism when compared to other baryogenesis scenarios. In the context of 
supersymmetric extensions of the SM it is possible (although not always simple) to combine 
the study of leptogenesis and neutrino physics with other physical phenomena like flavor- 
violating decays [18, 32], In general, this analysis does not give us definite answers, yet it 
may help to discriminate among certain neutrino mass models. 

Recently, a considerable amount of work has been done aiming at relating viable lepto- 
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genesis to all the available low-energy neutrino data coming from solar, atmospheric and 
reactor experiments [33]. Roughly speaking, two different approaches to the problem are to 
be found. The first one is based on the computation of the baryon asymmetry as a function 
of the lightest heavy Majorana neutrino mass M\ , the CP asymmetry eq and the so-called 
effective neutrino mass rhi = (m' D m D )u/Mi [34], By solving the Boltzmann equations, this 
kind of analysis provides valuable information on the ranges of these parameters that lead to 
an acceptable value of Y B . The weak point of this procedure lies on the fact that the input 
information depends on quantities which are not sensitive to the full structure of Y v and Mr 
and, therefore, no further conclusions can be drawn about the class of models which can 
lead to acceptable values of the input parameters referred above. In fact, the values of Mi, 
m i and eq should not be taken as independent parameters. The second approach is based 
upon initial assumptions on the structure of Y v and Mr at high energies which are fixed 
by recurring to theoretical arguments like for example grand unified theories or flavor sym¬ 
metries. Although in this framework some generality is lost, it has the advantage that one 
can compute the generated baryon asymmetry and, simultaneously, perform a low-energy 
neutrino data analysis. It is precisely the latter approach that we shall follow in the present 
work. 

In this section we present a class of minimal scenarios for lcptogenesis based on the trian¬ 
gular decomposition of Y v given in Eqs. (32) and (33). Namely, we would like to answer the 
question of how simple can the structure of the Dirac neutrino Yukawa coupling matrix be in 
order not only to get an acceptable value of the baryon asymmetry but also to accommodate 
the neutrino data provided by the atmospheric, solar and reactor neutrino experiments. In 
particular, we require the non-vanishing of the CP asymmetry generated in the decays of the 
lightest heavy Majorana neutrino, since the final value of the baryon asymmetry crucially 
depends on this quantity. Throughout our analysis we shall also consider the predictions on 
the CP-violating effects at low energies. 

In the previous section we have seen that the Dirac neutrino Yukawa coupling matrix Y v 
can be decomposed into the product of a unitary matrix U and a lower-triangular matrix 
Ya (cf. Eq. 32). It was also shown that the CP asymmetries £j generated in the heavy 
Majorana neutrino decays do not depend on the matrix U. In the special case U — U it 
is possible to establish the connection between leptogenesis, low-energy CP violation and 
neutrino mixing, since the same phases affect these phenomena. We classify this scenario 
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as a minimal scenario for leptogenesis and CP violation in the sense that the CP-violating 
sources that do not contribute to leptogenesis are neglected. On the other hand, if U ^ t 
this connection is not trivial. Therefore, from now on we will consider the case U = 1 which 
implies the following simple structure for the Dirac neutrino mass matrix: 


nio — uYa = v 


V2i e* 021 
v 2/31 e' 031 


0 0 ^ 

2/22 0 

!/::2 S/33 y 


(37) 


Then, from Eq. (14) the CP asymmetry generated in the decay of the heavy Majorana 
neutrino Nj is 


£ j = 


1 

8vr (Ha)^ 


E Im K ff A)?,] /«. 


(38) 


where 


Ha = YlY A , = 

with f(x) dehned in Eq. (15). 

From Eqs. (37) and (39) we readily obtain 


(39) 


Im[(iP A )^] = y\ x y\ 2 sin(20 2 i) + 2 y 2 \ y - 2 2 2/31 2/32 sin 0 1 + y\ x y| 2 sin 6 2 , 

Im[(i 7 A ) 3 i] = ;4 2/33 sin (2 0 31 ), 

Im[(^A)L] = 2/32 2/33 sin ( 2 ^32) , (40) 


with 0 , = 02 i + 03 i - 032 and 6* 2 = 2 ( 0 31 - 0 32 ). 

All the information about neutrino masses and mixing is fully contained in the effective 
neutrino mass matrix M. u which is determined through the seesaw formula given by Eq. (7). 
In this case 


M v = 


v 

Mi 


/ 


y 2 n 


2 /n 2 / 2 ie 


i </>2i 


3*031 


0 i 021 


2/n2/3ie 

2/n2/2ie‘^ 1 2/lie 2 * 021 + yi 2 jf-, 2/2i2/3ie* (031+</,2l) + 2 / 222 / 32 H e^ 32 

2/ii2/3ie^ 31 |/2i2/3ie i( ^ 31+02l) + 2 / 222 / 32 ^ e^ 32 y^e 2 ^ 31 + y 3 2 3 f| + y 3 2 2 f| e 2 ^ 32 




) 


(41) 


It follows from Eqs. (38)-(41) that, in principle, one can obtain simultaneously viable values 
for the CP asymmetries Ej and a phenomenologically acceptable effective neutrino mass 
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matrix in order to reproduce the solar, atmospheric and reactor neutrino data. This can 
be achieved by consistently choosing the values of the free parameters y tl , M* and faj . 
Yet, a closer look at Eqs. (38)-(40) shows that there are terms contributing to £j which 
vanish independently from the others. This means that a non-vanishing value of £j can be 
guaranteed even for simpler structures for Y„, which can be obtained from Ya assuming 
additional zero entries in the lower triangle'. The results are given in Table 1 where we 
present the textures constructed from Ya by neglecting one (textures 1-111) and two (textures 
IV-VI) off-diagonal entries. The form of the effective neutrino mass matrix AA V and the 
expressions for the CP asymmetries for each case are also given 7 8 . 

Let us first discuss textures IV-VI. For these three textures the effective neutrino mass 
matrix JA V predicts a complete decoupling of one light neutrino from the other two. This is 
in disagreement with the available neutrino data which indicates that, in the framework of 
the LMA solution, only one neutrino mixing angle should be small, namely 9 \ 3l instead of 
two as predicted by textures IV-VI. Furthermore, texture VI predicts a vanishing value for 
the CP asymmetry in the decay of the lightest heavy Majorana neutrino, implying Yb — 0 
in the case of hierarchical heavy Majorana neutrinos. Texture III can also be excluded on 
the grounds that it cannot account for the large solar angle. To illustrate this, let us write 
for this texture the matrix L given in Eq. (24) as 


L 


i rri D D r 


1/2 


1 z x 0 0 ^ 

^2 V 2 0 

\Z3 0 X 3 j 


(42) 


where zi, y 2 and x 3 are real and positive [26]. We also take z 2 and z 3 real to simplify the 
discussion. Considering first the hierarchical case Am@ ~ m 3 , a large atmospheric angle 
and m 2 < m 3 naturally arise if z 2 ~ z 3 > y 2 ,x 3 . More precisely, tan d 23 — z 2 /z 3 . In 
addition, this implies 

tan 0 13 ~ ^ = 1 13 , (43) 

V-2 + 4 


7 Notice however that the vanishing of diagonal elements in Ya would imply det (mu) = 0 and consequently, 
det (M. v ) = 0, leading to the existence of massless light neutrinos. 

8 In commonly used language, textures I-III and IV-VI belong to the classes of four and five texture zero 
matrices, respectively. For a complete discussion on seesaw realization of texture-zero mass matrices see 
e.g. [35] and for its implications at low energies see e.g. [36]. 
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TABLE I: Minimal textures based on the triangular decomposition of mp and their respective light neutrino mass matrix and CP- 
asymmetries £\ and £ 2 . 


Texture 


I 


II 


III 


IV 


V 


VI 


V A 


M * 


3/11 0 0 

3/21 e 1 ^ 21 3/22 0 

\ 0 2/32 e l ^ 32 2/33 y 


3/11 

0 


0 0 

3/22 0 


y 3/31 e* ^ 31 1/32 e l ^ 32 3/33 j 


( 3/11 0 0 ^ 

3/21 e 1 ^ 21 1/22 0 

y 3/31 e 1 ^ 31 0 2/33 y 


/ 3/11 o o \ 

3/21 e 1 ^ 21 3/22 0 

y 0 0 3/33 y 


/ 2/11 0 0 ^ 

0 2/22 0 

y 2/31 e 1 ^ 31 0 3/33 y 


/ 


2/ll 


2/112/21 e 


* 021 0 \ 
2/113/21 e 4 ^ 2 i 2/§j e 2i ^ 21 + ^ 3/| 2 2/22 3/32 ^ e* ^ 32 

0 1799 1799 Ml e i <?f>32 17oo Ml 4- 17 ? Ml ft 2 * ^ 32 , 

u y 22 y 32 m 2 e ^33 m 3 + ^32 m 2 e / 


2/11 


1,2 Ml 

*22 m 2 


3 /n 3/31 e l ^ 31 
3/323/22 Hje 1 ^ 32 


^ 3/113/31 e^ 31 2/32 3/22 ^ e^ 32 2/fi e 2 ^ 31 + 2/| 3 ^ + 2& ^ e 2 i <*> 32 J 


( 3 / 2 j 2/112/21 e *^ 21 2/112/31 e *^ 31 ^ 

2/112/21 e 1 ^ 21 2/|i e 2 *^ 21 + 3/| 2 ^ 3 / 21 3/31 e l ^ 21+</>31 ) 

^ 2/112/31 e 1 ^ 31 2/212/31 e l ^ 21 +^ 31 ) 2/§ x e 2 * ^ 31 + 2/33 ^ y 


/ 


2/ii 


2/112/21 e 


i 021 


2/112/21 e 1 ^ 21 2/|i e 21 * 2 i + 2/| 2 ^ 0 


V 


„2 Mi , 

i /33 M 3 / 


2/ll 


„2 Ml 
"22 m 2 


2/112/31 e 

0 


* <A31 \ 


\ 2/112/31 e^ 31 0 2/31 e 2; ^ 31 + 2/33 ^ / 


Stt (2 /| 2 +2/32 ) 

= 2/312/32 sin[ 2(</> 3 2 - 031 )] /21 - 2/312/33 sin( 20 3 i)/ 3 i 
®/i( 2 /ii + S/31) 

= 2/312/32 sin[2(0 3 i - </» 32 )] fi2 - 2/322/33 sin( 2^ 32 )/32 

87r (2/i 2 + 2 / 32 ) 

j/212/22 sin(2<^2i)/2l + 2/312/33 sin(20 3 i)/ 3 i 


£l = — 


871(2/1! +2/21 +2/31) 


. 2/212/22 sin ( 2 </> 2 l) t 

o 2 /12 

8712/22 

j/21 2(22 sin(2<^2l) 
811 (2/n + 2/21 ) 

. 2/212/22 sin(2 021) , 

o 2 /l2 

8tt y| 2 


-V . 
-I 


2/31 2/§3 sin(2(/>3i) 

“ MI&+I&) /31 


£9 = 0 


( 2/11 0 0 ^ 
0 1/22 0 
y 0 2/32 e* ^ 32 1/33 y 


^ 2/n 


0 2/22 Mj 


2/22 2/32 ^e^ 32 


y 0 2/22 2/32 ^ e i ^ 32 2/33 ^ + 2/32 ^ e2 * / 


£ VI = 0 


£9 = - 


2/32 2/33 sin (2 </> 32 ) 
Sir (2/22 + 2/32 ) 


/ 3 2 



which has to satisfy the CHOOZ bound, £13 < 0.2 [37]. The solar angle is approximately 
given by 


tan( 20 i 2 ) ~ 2 ti 3 


yl + %1 


< 


2 ti3 (3 


= B V 2 (f3, tiz ), 


(44) 


f 2 yl + f 3 xl-tu(f 3 yl + f 2 xl) i-ti 3 /? 2 
where (3 = max(tan 0 2 3 , cot 6 * 23 ). The upper bound -B 12 (/?, ^13 = 0.2) is an increasing function 
of f3. Experimentally 0.7 < (3 < 1.7, so that = E?i 2 (1.7, 0.2) = 4. This corresponds to 
tan 2 9 12 < 0.3, while the fitted LMA solution requires tan 2 0 V2 > 0.3. We thus conclude that 
type III cannot account for the observed large solar angle in the hierarchical case. Moreover, 
it turns out that it cannot also account for the pattern of M. v required by inverse hierarchy 
since, in this case, the atmospheric oscillation fit would require z\z 2 ~ Z 1 Z 3 ~ y/ A m@ and 
z i-> z 2 i z 3 ^ v/A^|- Fi na Hy) the degenerate case can be accommodated only by tuning the 
elements of L. In view of this, from now on we will focus our analysis only on textures I 
and II. 

Another interesting fact which comes out from the observation of Table I is that the 
phase content for textures I-II can be further reduced. Indeed, it is straightforward to show 
that with only one phase in Ya it is possible to obtain a non-vanishing CP asymmetry eq. 

As it has been stated in Section II, the strength of CP violation at low energies is sensitive 
to the value of the CP invariant Jcp given by Eq. (12). For the textures I and II we get: 


Jc 


CP 


yhyhyhyhv 12 
MfMf L 

M 2 


(y ‘21 y .32 + y'n y 2 22 + y'u yL) sin( 20 21 ) - y \ 2 y\^- sin( 2 . 032 


2/33(2/11 + 2/2i)-TT sin [2 (<£21 - <£32)] 


11 


3 CP 


M, 

yh yh 2/312/32 y 12 
AmljAm^Aml 


.(2/22 2/31 + vh y 2 22 + y 2 n y 2 32 ) sin t 2 (032 - M] 


32 


+ ^22^33^ Sin(2d 32 ) - y 2 u y 2 3 ^ sin(20 31 ) 


(45) 


From Table I and Eqs. (45) it is clear that we can have both e 1 ^ 0 and Jcp ^ 0 with a 
single non-vanishing phase. 


On the relative sign between Yg and Jcp 

It has recently been pointed out in [24] that the relative sign between CP violation and 
the baryon asymmetry can be predicted in a specific class of models. In our framework, it is 
clear from Table I and Eqs. (45) that the relative sign between the low-energy invariant Jcp 
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TABLE II: Relative sign between the CP asymmetries £12 and the low-energy CP invariant Jcp for 
textures I and II considering the different heavy Majorana mass regimes. We assume M\, M 2 <C M 3 
and (j >32 = 0. The parameter x — 0.4 corresponds to the zero of the loop function f(x ) defined in 


Eq. (15). 



Mi 

77 < X 

M 2 

Mi 

x < —— < 1 

M 2 

m 2 

—— < x 

Mi 

m 2 

x < —- < 1 

Mi 

sgn (e\ ■ j} Jp ) 

+ 

+ 

+ 

- 

sgn (4 ■ 3} JP ) 

- 

+ 

- 

- 

sgn (e\ l ■ 3}} P ) 

- 

- 

- 

+ 

sgn ( 4 1 • Jj! p ) 

+ 

- 

+ 

+ 


and the CP asymmetries £j cannot be predicted without specifying the values of the heavy 
neutrino masses Mi and the Dirac Yukawa couplings 9 . This is mainly due to the fact that at 
least one of these quantities depends on both phases appearing in the Dirac neutrino Yukawa 
coupling matrix for textures I and II. Therefore, in order to establish a direct connection 
between the sign of the CP asymmetries £j and the low-energy CP invariant Jcp further 
assumptions are required. For instance, considering the case ^32 = 0 and Mi,M 2 -C M 3 so 
that the terms proportional to /31 can be safely neglected in £ 1 ^ (see Table I), we obtain 
from Eqs. (45) and Table I the relative signs given in Table II. 

It is well known that in the case of hierarchical heavy neutrinos, the sign of Y B is fixed 
by the sign of the CP asymmetry generated in the decay of the lightest heavy neutrino, let 
us say £l . For Y B to be positive, as required by the observations, one must have £ P < 0. 
Thus, the measurement of the sign of CP violation at low-energy could in principle exclude 
some of the heavy Majorana neutrino mass regimes presented in Table II. As an example, 
let us assume that Jcp < 0 . In this case, the requirement £p < 0 would exclude the mass 
regime M 2 /Mi < x for texture I and Mi/M 2 < x for texture II. 

Finally, it is interesting to note that the Dirac neutrino Yukawa matrices considered in 
[24] correspond to our textures I and II with y 33 = 0 and have the remarkable feature 
that the number of real parameters equals the number of masses and mixing angles. In 
this limit the heavy Majorana neutrino IV 3 completely decouples, rendering this situation 

9 For the sake of more generality, we do not assume the mass-ordering M\ < M 2 in this discussion. 
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phenomenologically equivalent to the two heavy neutrino case considered in [24], Namely, it 
can be easily seen that for y 33 = 0 the phase (f) 32 can be rephased away and consequently the 
CP asymmetries e lj2 and the low-energy CP invariant Jcp will only depend on the phases 
02 i and (j )31 for textures I and II, respectively. This means that the examples considered in 
[24] are special cases of our framework which is motivated by the condition that rrifj does 
not contain any phase that would only contribute to low-energy CP violation. It is clear 
from Table II that, in general, even by keeping only one phase the sign of £j ■ Jcp depends 
on the particular Yukawa texture and hierarchy between the heavy Majorana neutrinos. 


VI. EXAMPLES 

In this section we present some examples of the minimal textures discussed in Section V 
and proceed to the study of their implications to low-energy physics as well as to the com¬ 
putation of the baryon asymmetry Yr through the numerical solution of the Boltzmann 
equations as described in the Appendix. We will only consider cases that lead to the LMA 
solution of the solar neutrino problem, which means that the neutrino mixing angles and 
the squared mass differences lie in the typical ranges: 

2.5 x 10" 5 eV 2 < A m% < 3.4 x 10' 4 eV 2 , 0.24 < tan 2 0 Q < 0.89 , 

1.4 x 10" 3 eV 2 < Am@ < 6.0 x 10” 3 eV 2 , 0.40 < tan 2 9@ < 3.0 , (46) 

with the solar and atmospheric mixing angles being identified as 9 Q = 9 1 2 , 9@ = 9 2 3, respec¬ 
tively, in the standard parametrization [13]. The 9 13 mixing angle is at present constrained 
by reactor neutrino experiments to satisfy | sin 9\ 3 \ = \U e3 \ < 0.2 [37]. 


A. Hierarchical heavy Majorana neutrinos (M\ < M 2 <C -M3) 


The first example is a realization of the texture I given in Table I with cp 32 = 0. The 
entries of the Dirac neutrino Yukawa coupling and the right-handed neutrino mass matrices 
are chosen to be of order: 


Yu = 


Ad 

v 


g 6 21 


0 0 
e 4 0 
e 4 1 


\ 


, Dr — Ar diag (e ,e , 1 ), 


(47) 
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where e < 1 is a small parameter. For our numerical estimates we consider the typical Dirac 
neutrino mass scale to be ~ 100 GeV, which corresponds approximately to the top quark 
mass at the GUT scale [38]. The neutrino mass matrix M. u is then given by: 


M v 


e 2 A r 


( e 2 
e 2 e*^ 21 

0 


e 2 e^ 21 
1 + e 2 e 2i<t>21 

1 


0 \ 


(48) 


In this particular case, considering for the moment 0 2 i = 0, one gets 10 : 


mi ~ 


A 2 , 

-P- (2 

2A fi l 




A 2 


m 2 — 


D 


2A 


R 


(2 + V2) 


m " 7T R (2 + e2) ' 


(49) 


leading to 


A — 


2v/2A 4 


° , Am. 2 2 ~ 4Ad 


32 ' 4 A| 


(50) 




The requirement 


AttTq 


Ami 


Am 2 j 


LMA 


Am| 2 C2 ' 


(51) 


forces e to be in the range 0.3 < e < 0.7. For the neutrino mixing angles we have: 


tan 2 0 i 2 — 1 


2^2 


, tan 2 0 23 ~1 , \U e3 \ ~ 


2 \/2 ’ 


(52) 


which are in good agreement with the data taking into account the range of e. The matrix 
Vl which corresponds to the left-handed rotation involved in the diagonalization of mo can 
be viewed as the equivalent in the lcptonic sector of the quark mixing matrix Vqkm- This 
matrix is obtained by diagonalizing the Hermitian matrix Y v Yj, 


Vr ~ 


1 - e 8 /2 
-e 4 + e 12 /2 


e 4 - e 12 /2 
1 - e 8 /2 


A 2 


—e L 


(53) 


This is in fact an interesting result in the sense that we are in the presence of a typical large 
mixing out of small mixing situation, where the large neutrino mixing is generated through 


10 We have checked that the light neutrino masses m, and mixing angles 9, j are not sensitive to the phase 
021 • 
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the seesaw mechanism [29, 39]. The scale Ar, or equivalently the mass M 3 of the heaviest 
Majorana neutrino, is determined by requiring Am| 2 to be in the experimental range given 
in (46). We find 


M 3 = A R ~ 


2 A l 


2 5 / 4 A 


D 


~ 8 x 10 15 GeV, 




(54) 


o 


and Eq. (47) implies Mi ~ 1.3 x 10 11 GeV and M 2 ~ 8.4 x 10 12 GeV, for e = 0.4. 
From the expression of the CP asymmetry e\ given in Table I, we obtain 


3 yl, l)l 2 sin (2 02i) Mi 

£ 1 =-~ 


3 A 2 c 10 

D sin(2 <t> 21 ) ~ 10 -6 sin(2 0 2 i) , 


(55) 


167T ( y\ x + y\ x ) M 2 327t v 2 

which is maximized for 0 2 1 = 7t/ 4. The value of the CP invariant Jcp can be computed 
from the expression 

M 2 


Jcp 


yh vii vh vh v 12 

MfM|Am|iAm|iAm| 2 |_ 


I/211/32 + 2 /n(2/22 + 2/32) + 2/33(2/11 + y 2 i)^ 


sin (20 


21, 


(56) 


which is obtained from Eq. (45) setting 0 32 

3 e 2 sin (2021) 


Jcp - 


32 


= 0. The result is straightforward in this case 
~ 1.1 x lO~ 2 sin(2 02i), (57) 


rendering visible CP-violating effects in the next generation of neutrino factories for 0 2 1 — 
7t/ 4. It is interesting to note from Eqs. (55) and (57) that the dependence of £1 and Jcp on 
the phase 0 2 i is such that both quantities are simultaneously maximized. Notice also that 
£1 and Jcp have the same sign, as it should be according to Table II. 

All the results presented above are confirmed by the full numerical computation presented 
in Fig. 1, where we have randomly included 0(1) coefficients 11 for the non-vanishing entries 
of Y„ and taken 0 2 i = 7 t/ 4. The spreading of the points in the figure, due to the random 
variation of the coefficients, shows that the textures are stable under these perturbations. 

In order to compute the value of the baryon asymmetry we proceed to the numerical solution 
of the Boltzmann equations as described in the Appendix, taking the initial conditions: 
Y£ = 0, Yg_ L = 0. The results are presented in Fig. 2 where we have plotted Y Nl , Y N2 and 
Yb as functions of the parameter z = M\/T for given values of the CP asymmetries, heavy 


11 For illustration, we have taken the 0(1) coefficients in the range [0.9,1.3]. We notice however that the 
results are not too sensitive to slight variations of this range. 
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FIG. 1: Results of the numerical computation for the textures presented in Eq. (47) for Y v and 
Mr with the CP-violating phase </> 2 i = 7 t/4 . The heavy Major ana neutrino masses M t and the CP 
asymmetries £\. e 2 are plotted as functions of the texture parameter e. The dependence of \U e %\ 
and Jcp on e is also shown. The dotted areas satisfy the neutrino constraints (46) for the LMA 
solar solution. 


Majorana neutrino masses and the lightest neutrino mass. The predicted value for the final 
baryon asymmetry is Yr — 9 x ICR 11 , which is inside the observational range (1). 

Our next example is a particular case of the texture II presented in Table I with </> 32 = 0. 
The Dirac neutrino Yukawa coupling and heavy neutrino mass matrices are chosen in this 
case to be of order 


Y — —— 
1 


A D 

V 


e 

0 

^5 <^31 


\ 


, D r = A r diag (e 10 , e 8 ,1), 


(58) 
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FIG. 2: Results for the distributions Y/v 15 Y/v 2 and for the baryon asymmetry Y R as functions of 
z = M\/T, obtained from the numerical solution of the Boltzmann equations. The plot refers to 
the texture for Y u considered in Eq. (47) with cf >21 = vr/4 and e = 0.5. The values of the heavy 
Majorana neutrino masses Mi, the CP asymmetries eq 2 and the lightest neutrino mass m\ are 
consistent with the ones presented in Fig. 1 for this value of e. The value predicted for the final 
baryon asymmetry is Yb — 9 x 10 -11 . 


leading to the following light neutrino mass matrix 


M v 


A 2 

e 2 A R 


( 


0 

^2 31 


0 e 2 e irt>31 

1 1 
1 1 + e 2 (l + e 2 *^ 31 ) y 


(59) 


The predictions for the Arne’s and neutrino mixing angles are similar to the ones given 
in Eqs. (50) and (52). Moreover, from Eqs. (54) and (58) and for e = 0.4, we get the 
following heavy Majorana neutrino masses: Mi ~ 9 x 10 11 GeV, M 2 ~ 6 x 10 12 GeV, M 3 ~ 
8 x 10 15 GeV. The left-handed matrix Vl is given in this case by: 


Vl — 


-e 

no 


10 


—e 

1 

-e 6 


10 


no 


\ 


(60) 


24 













FIG. 3: The same as in Fig. 1 for the textures considered in Eq. (58) and 0 31 = 7 r/ 4 . 


Finally, the CP asymmetry can be obtained from the expression of e * 1 presented in Table I 
with 032 = 0. Taking into account the form of Y v and Dp as in Eqs. (58), we obtain 
3 2 / 3 ! sin(20 3 i) / , Mi 0 M ] \ 3A 2 n e 8 


£1 = 


167 rivh + yh) 


2 Ml | 2 Ml \ a 1V D < 

V:i2 M 2 +V:i:i M :i ) - 32 7ru 2 


sin(2 03 1 ) ~ 6.5 x 10 b sin(2 0 3 i), 


(61) 


where the last estimate has been obtained assuming e = 0.4. It can also be shown that the 
CP invariant Jcp reads (see Eq. (45)), 

M 2 


Jcp = 


vh vi 2 y'ii yh v 12 


MfM 2 Am| 

which leads to the approximate result 

3 e 2 sin (2031 


Z 2 1 2(2 1 z\i 2 2 " 

y 221/31 + y n J22 +1/32) + Vh 2/33 


Jcp — 


32 \pl 


~ — 1.1 x 10 2 sin (2 0 31 ). 


sin ( 2031 ), (62) 


(63) 


Notice that £1 and Jcp have opposite signs in this case, once again in agreement with 
Table II. 

In Fig. 3 we present the same numerical analysis as in Fig. 1, but for the case where Y u 
and Dp are defined through Eqs. (58). We find good agreement between our approximate 
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z=M 1 /T 


FIG. 4: The same as in Fig. 2 for the example considered in Eq. (58) with foi = 7t/4 and e = 0.4. 
The values of the heavy Majorana neutrino masses Mj, the CP asymmetries £ 1.2 and the lightest 
neutrino mass mi are consistent with the ones presented in Fig. 3 for this value of e. The predicted 
value for the final baryon asymmetry is Yb ~ 8 X 10~ n . 

analytical results and the numerical ones. The integration of the Boltzmann equations is 
plotted in Fig. 4. The value for the final baryon asymmetry is Yb — 8 x 10 -11 . 


B. Two-fold quasi-degeneracy (Mi ~ M2 M3) 

As an example, let us consider the texture of type I given in Table I and assume that 
Mi ~ M2 -C M3. The Hermitian matrix H A = Y^Y A is given by 


( 


H a = 


2/11 + 2/21 2/212/22 e + 21 


0 


V 


2/212/22 e 

0 


2/22 + 2/32 2/32 2/33 e" 


2/32 2/33 e 


*032 


2/33 


\ 


/ 


(64) 


Taking into account the requirement of Eq. (20) we can get the range of validity of the 
parameter 5n- 


> ^ max {2/11 + 2/21 , 2/22 + 2/32} 


(65) 
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In the limit Sn —> 0, Eq. (48) is recovered. Since Sn -C 1 the results for the neutrino 
masses and mixing parameters are practically the same as the ones given in Eqs. (49)-(52). 
The same is expected for the estimate of the CP invariant Jcp defined in Eq. ( 66 ). The 
main differences reside obviously in the heavy Majorana neutrino mass spectrum and on the 
values of the CP asymmetries, since now the relation Mi <C M 2 is no longer valid. From 
Eq. (23) we obtain 


e 18 A 2 d sin(2 0 2 i) 1.3x10 12 . . , 

£l ~ -55—1-5---7-sin(2 02 i) , 

32 7 r On v z On 

e 20 A 2 D sin (2 0 21 ) 1.1 x 10" 13 . . 

£2 -- 55—7 - 5 --- 7 -sm (2 0 21 ) • 

32 7 T On On 


The results of the numerical integration of the Boltzmann equations for this case are pre¬ 
sented in Fig. 5. A realistic value for Y R is also generated in this case. 

Before we end this section, it is worthwhile to comment on the possible effects of quantum 
corrections to the effective neutrino mass matrix AA V . This discussion turns out to be 
relevant since in the examples considered above we have taken the effective neutrino mass 
matrix at the heavy neutrino decoupling scale. Although a detailed treatment would require 
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z=M 1 /T 


FIG. 5: The same as Figs. 2 and 4 for the quasi-degenerate case M\ ~ M -2 considered in Eq. (67) 
with (f> 2 i = 7 t/4 , e = 0.3 and <5 jv = 5 x 10“'. The predicted value of the baryon asymmetry is 
Y b ~ 10 -10 . 

a renormalization group analysis for the effective neutrino mass operator, one can employ the 
simple analytical treatment considered by many authors in the literature (see for example 
Refs. [22, 40]). Following this, we recall that the effective neutrino mass matrix at the 
electroweak scale rriz can be related to the one at the heavy neutrino decoupling scale A B 
as 


M. imz) oc diag(l, 1,1 + e T )M u (^r) diag(l, 1,1 + e T ), 


(71) 


where, neglecting the running of the charged lepton Yukawa couplings, 

327r 2 \rriz 


(72) 


Taking y T = m T /v and Ar < 10 16 GeV, we obtain e T < 10 5 . Considering M. v (Ar) as given 
in Eq. (48) we would get from Eq. (71): 


(1 -O' 
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, tan 2 9 12 — 1 - -y= — V2 e T , 
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2V2 \ 1 _ 2 & 


(73) 












Since e r <C 1, the results given by Eqs. (51) and (52) are not significantly altered. The same 
conclusions are drawn for the example considered in Eq. (59). Thus, we conclude that the 
effects of quantum corrections due to the renormalization of the neutrino mass operator can 
be safely neglected in our case. 

VII. CONCLUDING REMARKS 

We have analyzed, in the context of the minimal seesaw mechanism, the link between 
leptogenesis and CP violation at low energies. In particular, it was shown that, in order to 
present a thorough discussion on this question, it is convenient to work in the WB where 
both the charged lepton and right-handed Majorana neutrino mass matrices are diagonal 
and real, and to write, without loss of generality, the Dirac neutrino mass matrix as the 
product of a unitary matrix and a lower triangular matrix. From the analysis of the phases 
that contribute to leptogenesis and low-energy CP violation, we have identified a necessary 
condition which is required in order to establish a link between these two phenomena. We 
have studied a class of models which satisfy the above necessary condition in the simplest 
way, namely those where the Dirac neutrino mass matrix is of the triangular form. By 
choosing this structure the number of physical parameter in the theory is reduced then 
enhancing its predictability. In this case there are only three CP-violating phases which 
contribute both to leptogenesis and CP violation at low energies. We have then studied the 
minimal scenarios where a correct value of the baryon-to-entropy ratio can be generated, 
while accounting for all the low-energy neutrino data in the context of the LMA solution. 
Moreover, the examples considered in Section VI predict the existence of low-energy CP- 
violating effects within the range of sensitivity of the future long-baseline neutrino oscillation 
experiments. In fact, it is a remarkable feature of these scenarios that the solutions viable for 
leptogenesis are precisely those which maximize Jcp ■ The question of relating the observed 
sign of the baryon asymmetry to the sign of the leptonic CP violation, measurable at low 
energies through neutrino oscillations, was also considered. Namely, we have concluded 
that, within the minimal scenarios presented, this relation crucially depends on the heavy 
Majorana neutrino mass spectrum. We remark that a full discussion of this aspect requires 
the computation of the BAU since, besides the prediction of the relative sign between the 
BAU and Jcp-, the determination of Yp is of extreme importance to infer about the viability 
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of a given model. 
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APPENDIX: BOLTZMANN EQUATIONS 

The computation of the cosmological baryon asymmetry Yb involves the solution of the 
full set of Boltzmann equations which governs the time evolution of the right-handed neutrino 
number densities Y N . and the generated lepton asymmetry Y L . These quantities depend not 
only on the physics occurring in the thermal bath but also on the universe expansion. In the 
SM framework extended with heavy right-handed neutrinos the physical processes relevant 
to the generation of the baryon asymmetry are typically the N t decays and inverse decays 
into Higgs bosons and leptons, the A L = 1 scatterings involving the top quark and A L = 2 
scatterings with virtual N t and A/ — Nj scatterings [41, 42], The initially produced lepton 
asymmetry Y L is converted into a net baryon asymmetry Yb through the (B + /./-violating 
sphaleron processes. One finds the relation 

Y b = ZY b _ l = ^jY l . (A.l) 

For the SM with three heavy Majorana neutrinos £ ~ 1/3 and therefore Y B — Y B - l /3 — 
- Y l /2 . 

The Boltzmann equations for the N t number densities and the Y B -l asymmetry in terms 
of the dimensionless parameter z = Mi/T, where T is the temperature, can be written in 


30 



the form: 




NN ( Y nYn 3 

ij l y eq Y eq 


dV B -L 

dz 


Hs(z) 


^-i) + 


W Y B-L 


(A. 2 ) 


where H is the Hubble parameter evaluated at z — 1 and s(z) is the entropy density given 
by 


1 47T 3 g* Mf 
45 Mp 


s(z = 


2n 2 g* Mf 
45 £ 3 ’ 


(A.3) 


respectively. Here, g* ~ 106.75 is the effective number of relativistic degrees of freedom and 
Mp ~ 1.2 x 10 19 GeV is the Planck mass. The equilibrium number density of a particle i 
with mass rn l is given by 


y- q (z) = 


45 g { (rrii\ 2 (rn % z 


4 Mg* \Mi 


Mi ) ’ 


(A.4) 


where g t denote the internal degrees of freedom of the corresponding particle (g^. = 2 , 
g( = 4) and K n (x) are the modified Bessel functions. The q’s are the reaction densities for 
the different processes. For the decays one has 


7 ? = 


MiMf 




(A.5) 


The reaction densities for the A L = 1 processes with the top quark and for the A L = 2 
scatterings N t — N } can be written in the following way: 


7? = A” + . %f' = 7$v, + 7A, ' 


(A. 6 ) 


respectively. Finally, 7 w accounts for the washout processes: 


^ = E ^ + §A>+27! 


« +2 7 « + 2 7 f. 


(A.7) 


Each of the above reaction densities can be computed through the corresponding reduced 
cross section a^\x): 


Mi) (z) Ml 

64 7 r 4 z 


(x) \fx Ki (z y/x) dx , 


(A. 8 ) 


(ml+mD/M* 
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where x = s/Mf , with s being the center-of-mass energy squared and m a ,b are the masses 
of the initial particle states. All the relevant reduced cross sections are summarized below. 
For a more detailed presentation the reader is addressed to Ref. [42], 

We write the reduced cross sections as functions of the parameter x, the Hermitian matrix 
H u = YjY„ and the quantities aj and Cj defined as 

a> (H„)% 


a j — 


Mi 

M 1 


c j = 


El 

M 1 


64 7 r 2 


where Tj is the decay rate defined in Eq. (19). 

For the A L = 1 processes involving interactions with quarks we have 
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where Oh (introduced to regularize the infrared divergencies) and a u are given by 
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(A. 10) 


(A. 11) 


47 T 47 TV 2 

respectively. The mass parameter // is chosen to be fi = 800 GeV as in Refs. [41, 42], 
The reduced cross sections relative to the N t — Nj scatterings are: 
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where the functions L l3 and X tJ read as: 


(A. 14) 


= 2 In | 
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(A.15) 


Since the heavy Majorana neutrinos are unstable particles, one has to replace the usual 
fermion propagators by the off-shell propagators Dj(x ) in the computation of the amplitudes 
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for the right-handed neutrino mediated processes: 


x — ( 1 ; 


Dj Dj(x) (x — dj) 2 3 + a,jCj 

For the A L = 2 processes the reduced cross sections read then 
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In our analysis we have computed numerically the reaction densities through Eq. (A.8) 
and the above definitions of the reduced cross sections. Nevertheless, useful analytical 
approximations can be obtained for specific temperature regimes [42], 
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